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Abstract: Coal mining has changed the hydrogeological conditions of river basins, and studying how the 
relationship among different types of water body has changed under the influence of coal mining is of 
great significance for understanding the regional hydrological cycle. We analyzed the temporal and spatial 
distribution of hydrochemical properties and environmental isotopes in the Hailiutu River Basin (HRB), 
China with a mixed model. The results showed that: (1) human activity (e.g., coal mining and agricultural 
production) causes considerable changes in the hydrochemical properties of surface water in and around 
the mining areas, and leads to significant increases in the concentrations of Nat and SOF ; (2) 
precipitation is the main source of water vapour in the HRB. The transformation between surface water 
and groundwater in the natural watershed is mainly affected by precipitation; and (3) in the mining areas, 
the average contribution rates of precipitation to the recharge of surface water and groundwater increased 
by 2.6%—-7.9% and 2.7%—9.9%, respectively. Groundwater in the Salawusu Formation constitutes up to 
61.3%-72.4% of mine water. Overall, this study is beneficial for quantifying the effects of coal mining on 
local hydrological cycles. The research results can provide a reference for local water resources 
management and ecological environment improvement. 
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1 Introduction 


The transformation among precipitation, surface water, and groundwater is important for 
understanding the hydrological cycle of arid or semi-arid regions. Coal mining in a river basin 
damages its aquifer structure, thereby complicating the transformation among precipitation, 
surface water, groundwater, and mine water, and obscures the characteristics, mode, and changes 
in these water transfers. 

The transformation from river water to groundwater in a region constitutes a crucial step in its 
hydrological cycle. The spatial and temporal distribution and evolution of natural water bodies 
can be understood by identifying their hydrochemical components (Chen et al., 2020; Jampani et 
al., 2020; Lei et al., 2021). The water vapour source, hydrochemical composition, and influencing 
factors of the basin are analyzed, and the temporal and spatial changes in hydrochemistry under 
the influence of different environments are analyzed (Sun and Chen, 2018; Nofal et al., 2019; 
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Wagh et al., 2020). The "tracer" effect of environmental isotopes can be used to further track the 
migration pathways of different water components, and interpret the transformation between 
surface water and groundwater (Hao et al., 2019; He et al., 2019; Du et al., 2021). Experts and 
scholars have conducted researches using hydrochemistry and isotope technology. A qualitative 
analysis of groundwater hydrochemical characteristics, pollution levels, and human health risk 
assessment that used temporal and spatial variations in groundwater hydrochemistry was 
conducted in western India (Kadam et al., 2021a, b). Hydrochemical and isotopic methods were 
used to study the mineralization processes of groundwater in the Walawe River Basin of Sri 
Lanka (Jayawardana, 2020). Mirzavand et al. (2020) analyzed the mechanism of groundwater 
salinization in the Kazan Basin of Iran with hydrochemical and isotopic methods. A study of 
groundwater in the Teboursouk Basin (Ayadi et al., 2018) in northwestern Tunisia showed that the 
ion exchange process, rock dissolution, and precipitation infiltration are the main factors that 
influence groundwater mineralization and the spatial distribution of groundwater chemistry. 
Using hydrochemical and isotopic methods, Rezaei et al. (2018) found that local precipitation is 
the main source of groundwater in the high-altitude area of Iran. Seasonal factors play a pivotal 
role in the transformation between surface water and groundwater in the arid and semi-arid 
regions of China (Wu et al., 2016; Chai et al., 2021; Song et al., 2021). Economic and social 
progress has recently heightened the impact of domestic and industrial activity on the 
transformation between surface water and groundwater. The aforementioned researches show that 
hydrochemistry and isotopes can be used to analyze both the changes in surface water and 
groundwater quality, and the supplement and discharge relationships among different types of 
water body. 

Coal is a major source of energy that propels social advancement. However, extensive 
excavation of coal from mines increases regional hydrological cycle. Several researchers have 
examined the effects of coal mining on the transformation between groundwater and surface 
water (Wang et al., 2016; Qu et al., 2018; Hao et al., 2019). However, the temporal and spatial 
variations in the transformation among precipitation, surface water, groundwater, and mine water 
in mined and unmined areas of river basins are not yet understood. Research has shown that 
large-scale excavation destroys original topographical and rock structures (Batsaikhan et al., 2017; 
Svobodova et al., 2019; Wang et al., 2020). For example, coal mining in eastern Donbas, Ukraine, 
has worsened surface settlement and rock deformation, resulting in a considerable decrease in the 
discharge of springs and surface water (Gavrishin, 2018). In addition, mining activity in the East 
Singhbhum District of the Jharkhand state in eastern India has led to surface settlement and the 
formation of cracks in several surface areas, which has altered groundwater discharge and 
lowered water levels in the shallow aquifers of the surrounding areas (Singh et al., 2018). Coal 
mining in the mining areas of northwestern China has brought about structural damage to the 
fractured and porous aquifers, and therefore caused a concentrated discharge of groundwater 
towards the goafs, resulting in increases in the recharge of groundwater by surface water (Guan et 
al., 2019; Guo et al., 2019; Han et al., 2020). 

By analyzing the spatial variation in hydrochemical properties and environmental isotopes, this 
study determines the changes of transformation among different types of water body and the 
influencing factors in natural watersheds and coal mining areas, and quantitatively identifies the 
amount of water transformed using the isotope mass balance equation. 


2 Study area and methods 


2.1 Study area 

Situated in the transitional zone between the Maowusu Desert and the Loess Plateau in northern 
Shaanxi, China, the Hailiutu River Basin (HRB) = (38°01'20"—38°50'42'N, 
108°38'21"—109°24'19"E) encompasses an area of 3326 km/?, and ranges in elevation from 1002 
to 1477 m a.s.l., with high elevation in northwest and low in southeast. The terrain consists 
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primarily of tidal flats and sandy lands, and its surface is covered mostly by aeolian sand, which 
is mainly composed of medium and fine grains, and is highly permeable, therefore favoring 
precipitation infiltration. Below the surface, the aquifer exists in the Salawusu Formation, which 
features a lithology mainly composed of medium and fine sand grains with a high water storage 
capacity. This aquifer serves as a major source of water for industrial and domestic activity within 
the region. Located within the semi-arid interior of northwestern China, the HRB has a temperate 
continental arid and semi-arid monsoonal climate characterized by high temperatures and 
precipitation in summer, and low temperatures and dry air in winter. Precipitation in this region 
considerably varies among seasons. 

According to the soil moisture characteristics, groundwater aquifers in the study area from top 
to bottom are as follows: Quaternary Holocene eolian and alluvial aquifers, upper Pleistocene 
Salawusu Formation pore aquifer, pore fissure phreatic water and confined water aquifer of the 
Lower Cretaceous Luohe formation, Aquifuge of the Anding formation of the Middle Jurassic 
system, the Middle Jurassic Zhiluo Formation confined aquifer group, and the Middle Jurassic 
Yan'an formation confined aquifer group. The loose pore aquifer of the Salawusu Formation is 
shallowly buried, and is an important ecological water source and water supply source. 


2.2 Sampling and testing 


The Hailiutu River Basin (HRB) is composed of the Hailiutu River (HR) and Yingdiliang River 
(YR; Fig. 1). Based on hydrogeological and field surveys, and the geographic locations of coal 
mines in the area, we divided the study area into two parts: the natural watershed, and mining area 
and downstream of the HRB. The natural watershed is located in the upstream of the HR and YR, 
and its population is relatively small. The conversion between different types of water body is less 
affected by human activities and coal mining in this region. The mining areas and downstream of 
the HRB mainly include Yingpanhao (K1), Dahaize (K2), and Balasu (K3 and K4) coal mines. 
Figure 1 shows the sampling sites. A total of 118 water samples were collected during the normal 
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Fig.1 Location of the sampling sites in the Hailiutu River Basin 
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water season (May), wet season (August), and dry season (October) of 2019. A total of 32 water 
samples were collected during normal water season (groundwater samples: H1—H5, H8, H12—H16, 
H18, Y1, and Y2; surface water samples: h1—h4, h6, h7, and yl—y3; mine water samples: K1—K4; 
and precipitation samples: JS1—JS4), and 86 water samples were collected during wet and dry 
seasons (groundwater samples: H1—-H18, Y1, and Y2; surface water samples: hl—h11 and yl—y4; 
mine water samples: K1—K4; and precipitation samples: JS1—JS4). Each water sample collected 
was placed in a 250-mL brown polyethylene sampling bottle that had been previously cleaned. 
During this process, the sampling bottle was rinsed three to five times. Each groundwater sample 
was bottled 1 to 3 min after extraction. Each sampling bottle was completely filled with the 
corresponding water sample, and no bubbles existed for preventing from evaporation. The 
sampling bottles were sealed with Parafilm, and placed in an insulated box at 0°C-4°C, after 
which they were transported to the laboratory for testing. 

An HI98129/H198130 portable multiparameter water quality meter (HANNA instruments, Italy) 
was used to measure the pH, degree of mineralization (i.e., content of total dissolved solids (TDS)), 
and electrical conductivity of the water samples. In addition, sampling site information was 
recorded. The concentrations of conventional ions (K+, Ca?*, Nat, Mg?*, CI, NO;, and SO; ) 
were measured using a Dionex Aquion Ion Chromatography (Thermo Fisher Scientific, USA). 
The concentration of HCO; was measured through hydrochloric acid titration with methyl orange. 
The analytical accuracy of cations and anions was determined by the ion balance error (IBE). Ions 
were calculated with the unit of mg/L, and the error of IBE observations was +5%. The 
concentrations of heavy metal elements (Cu, Mn, Al, Ni, Pb, Cu, Ge, and Ag) were determined 
using an inductively coupled plasma source mass spectrometer (Thermo Fisher Scientific, USA). 
Stable H and O isotopes (i.e., D and '8O, respectively) were measured with an LGR IWA-45EP 
liquid-water isotope analyser (IWA-45EP, Los Gatos Research Inc., USA) with testing accuracies 
better than +0.1% and +0.5% for the isotopic ratio of !8O (5!8O) to D (6D), respectively. The 
isotopic results are expressed in thousands of deviations relative to V-SMOW (Vienna standard 
mean ocean water), and the isotope ratio 6 is calculated as follows: 


Ream le Ry. 
5(%o) = sampl: V-SMOW , (1) 
Ry smow 
where Rsampte is the ratio of D/H or !8O/'6O in the water sample; and Ry-smow is the ratio of D/H or 


180/160 of V-SMOW standard water sample. 
2.3 Data analysis 


The hydrochemical types of the water samples were analyzed using a Piper trilinear diagram and 
Schukarev classification. A Piper trilinear diagram was drawn using AqQA software, and the 
Schukarev classification was calculated and analyzed in Excel 2013. 

The Gibbs diagram was used to analyze the factors that influenced the water hydrochemical 
types, and was drawn in Excel 2013. 

We established a multivariate linear mixed model according to the mass and concentration 
balance equation to quantitatively analyze the rates of contribution of atmospheric precipitation, 
surface water, groundwater, and mine water, which were calculated using IsoSource software. The 
formulae for these calculations are as follows: 

dm (%o)= fid1 + f2d2 + ... + frOns (2) 

fitfot...tfr=l, (3) 

where 6m is the 8D or 5!8O of the mixed water; 5, is the 5D or 8!8O of the n'* endmember; and f, 
is the mixing ratio of the n"" endmember. 


3 Results 


3.1 Hydrochemical characteristics 


According to the results of the main concentrations of surface water and groundwater in the HRB 
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(Table 1), mean TDS values of surface water in the study area were 470, 367, and 403 mg/L 
during the normal water, wet, and dry seasons, respectively, while the pH values of surface water 
did not significantly vary seasonally, and were generally weakly. In each season, Na* and HCO; 
were the dominant cation and anion in surface water, respectively. The maximum values of TDS 
and SO2.4 in each season were measured near the downstream of the Tuanjie Reservoir near the 
Yingpanhao mining area. This result indicated that drainage during coal mining period caused 
damage to the surface water quality of the HRB. The Na* content of surface water was higher in 
the normal water season, because agricultural production activity occurred from May to July in 
northern China, and the concentration of Na* in surface water increased due to the impact of 

agricultural production activities. 

Table 1 Major ions concentrations in surface water and groundwater of the Hailiutu River Basin 
Surface water 

Statistic 
Season" dex Cat Mg Nat Kt cl SO? HCO; NO; TDS 
(mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) (mg/L) 


Max 162.91 57.06 1762.11 14.15 318.44 2679.93 469.88 11.02 1258.00 9.07 


Min 4.70 078 29.97 0.35 33.21 98.90 181.23 3.49 212.00 8.24 
pia Average 46.74 19.25 788.37 5.70 133.14 907.40 274.81 7.44 576.00 8.43 
SD 49.98 19.25 582.12 4.90 98.59 949.17 94.39 280 381.00 0.28 
CV 1.07 1.00 0.74 0.86 0.74 1.05 034 0.38 0.66 0.03 
Max 150.43 70.24 464.34 6.71 69.52 609.52 276.42 12.25 997.00 9.04 
Min 60.70 27.08 17.34 0.84 10.50 41.10 198.93 0.05 176.00 8.01 
Wet average 99.37 41.57 134.58 265 31.13 171.53 245.69 425 42200 8.37 
SD 25.75 14.94 15433 2.16 17.35 164.46 24.08 3.92 277.00 0.27 
CV 0.26 0.36 1.15 0.82 0.56 0.96 0.10 0.92 0.66 0.03 
Max 262.82 76.73. 1100.39 18.27 134.36 1810.76 338.05 33.18 1090.00 8.83 
Min 62.85 4.41 54.22 1.76 17.01 48.87 22883 1.81 191.00 8.20 
Dry Average 131.77 39.98 374.63 7.77 57.12 506.39 275.14 841 473.00 8.53 
SD 56.25 28.59 338.33 4.94 36.01 570.91 29.65 881 291.00 0.15 
CV 0.43 0.72 0.90 0.64 0.63 1.13 O11 1.05 0.62 0.02 
Groundwater 
Max 122.86 20.13 600.50 4.59 107.12 174.63 268.22 226.08 258.00 8.16 
Min 2.86 0.47 10.25 034 14.26 19.55 144.98 1.61 93.00 7.55 
phi Average 26.80 7.76 82.39 1.09 32.68 74.78 185.29 51.94 169.00 7.86 
SD 32.19 649 15893 1.08 24.92 4310 34.09 61.61 4240 0.17 
CV 1.20 0.84 1.93 0.99 0.76 0.58 0.18 1.19 0.25 0.02 
Max 159.69 57.97 79.32 15.52 37.16 69.32 336.83 69.16 328.00 8.56 
Min 2.00 4.91 6.35 0.03 3.53 12.52 124.48 0.00 109.00 7.71 
Wet Average 68.03 25.01 26.00 1.31 11.18 25.46 201.60 16.09 176.00 8.06 
SD 35.22 14.35 19.68 3.46 9.07 1416 49.99 17.92 59.04 0.20 
CV 0.52 0.57 0.76 265 0.81 0.56 0.25 1.11 0.34 0.02 
Max 150.42 55.63 224.22 10.57 39.42 145.29 372.83 153.24 245.00 8.67 
Min 21.39 278 18.85 0.21 5.58 10.32 155.60 0.16 103.00 7.67 
Dry Average 78.38 21.42 59.85 212 13.11 43.36 229.40 30.53 168.00 8.23 
SD 32.77 15.85 50.91 2.69 8.43 35.53 58.16 44.45 47.28 0.26 
cy 0.42 0.74 0.85 1.27 0.64 0.82 0.25 1.46 0.28 0.03 


Note: Max, maximum; Min, minimum; SD, standard deviation; CV, coefficient of variation. The abbreviations are the same in Tables 3 
and 4. 
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TDS contents of groundwater were 173, 177, and 186 mg/L in the normal water, wet, and dry 
seasons, respectively, while the pH values varied within the weakly alkaline range. Evidently, the 
TDS content of groundwater was low, and did not significantly vary. Na* and Ca** were the 
dominant cations in groundwater, while HCO; was the dominant anion. The average HCO; 
concentration in groundwater was the largest, and the variation was relatively minor, which 
indicated that HCO; concentration of groundwater in the HRB accounted for the largest 
proportion, and was stable. Na*, Mg?*, CI, and SO; varied widely in the basin, which showed 
that the ion concentrations in different areas greatly varied, and was vulnerable to natural and 
human activities. 


3.2 Hydrochemical composition characteristics 


A Piper trilinear diagram was plotted to classify the hydrochemical types of surface water, 
groundwater, and mine water. In the unmined area (Fig. 2), HCO3-Ca-Mg was the dominant 
hydrochemical type of surface water, while HCO3-Ca-Mg were the main hydrochemical type of 
the groundwater. In comparison, in the mining areas and downstream of the TRB, the dominant 
hydrochemical types of surface water transitioned to HCO3-SO4-Na-Ca, while HCO3-Ca-Mg and 
HCO3-Ca-Mg were the main hydrochemical types of groundwater. There were minor differences 
in the hydrochemical types of the mine water, with HCO3-SO,4-Ca as its primary hydrochemical 
type. Analysis of the Piper trilinear diagram revealed similar values in surface water and 
groundwater sampling sites in the mining areas and downstream of the TRB, which can be 
ascribed to the higher population density and the presence of a large number of mines within this 
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Fig. 2 Piper trilinear diagram for different types of water body in the Hailiutu River Basin. (a), normal water 
season; (b), wet season; (c), dry season. Groundwater samples: H1-H5, H8, H12—H16, H18, Y1, and Y2; surface 
water samples: hl—h4, h6, h7, and yl—y3; mine water samples: K1—K4; and precipitation samples: J1—J4 for 
normal water season. Groundwater samples: H1-H18, Y1 and Y2; surface water samples: hl—-h11 and yl—y4; 
mine water samples: K1—K4; and precipitation samples: J1—J4 for wet and dry seasons. 
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area. Factors such as dewatering and drainage during mining as well as human activity lead to a 
stronger hydraulic connection between surface water and groundwater. The exchange between 
surface water and groundwater was more frequent during the wet season due to increased 
precipitation. As a consequence of dewatering and drainage operations during mining, SO; more 
strongly affected the hydrochemical property in the middle and upper reaches of the HRB, which 
were close to the Yingpanhao coal mine area. SO4-Na were the hydrochemical type of surface 
water (hl—h3) in this area, which can be primarily attributed to the abnormal increase SO; 

content caused by the prolonged residence time of the gypsum contained in surface water of the 
reservoir under construction in this area. 


3.3 Factors influencing hydrochemical composition 


The Gibbs diagram can facilitate analyses of the hydrochemical composition and mechanism of 
formation of natural watershed, and can directly reflect which of the principal 
factors—evaporation and concentration, rock weathering, or precipitation—govern their 
hydrochemical composition. In the Gibbs diagram, the logarithm of TDS (y-axis) was plotted 
against the Na*/(Na*t+Ca**) and Cl-/(CI-+HCO; ) ratios (x-axis). Figure 3 showed the Gibbs 
diagram for different types of water body in the study area. In the study area, rock weathering as 
well as evaporation and concentration collectively govern the hydrochemical types of surface 
water, while rock weathering primarily affected the hydrochemical types of groundwater and 
mine water. Compared with the unmined area, more surface water sampling sites in the mining 
areas and downstream of the HRB fall outside the Gibbs diagram, which suggested that human 
activity significantly impacted the hydrochemical types of surface water in this area. A greater 
number of surface water and groundwater sampling sites were found to deviate from the Gibbs 
model during the normal water season than during the other two seasons, which may be attributed 
to the greater impact of agricultural production activity on the hydrochemical types due to its 
higher frequency during the normal water season. 
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Fig. 3 Gibbs diagram for different types of water body. (al and a2), normal water season; (b1 and b2), wet 
season; (cl and c2), dry season. 
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The Gibbs diagram cannot determine the hydrochemical source and formation process of water. 
Therefore, the relationships between different ions can be determined by using the ion 
concentration of water samples, and the main ion sources of surface water, groundwater, and mine 
water in the study area can be analyzed. When (Ca*++Mg**)/((HCO,+SO; ) was greater than 1, the 
Ca** and Mg’ in water were mainly derived from carbonate dissolution; when this value was less 
than 1, the Ca** and Mg** in water were mainly derived from silicate dissolution. Most of the 
water samples in the ratio diagram between surface water and groundwater (Ca**+Mg**)/(HCO,;+ 
SO; ) in the natural watershed were distributed at the upper left of the 1:1 line, which showed that 
the main ions in waters of the HRB came from the dissolution of carbonate (Fig. 4). In the 
(Ca?*+Mg**)/SO; diagram, most of surface water and groundwater samples were distributed on 
the upper left side of the 1:1 line. The surface and groundwater samples were distributed on both 
sides of the 1:1 line in the normal water season, which indicated that Ca2*++Mg?* and SO; 
contents of the water in the normal water season were similar, while Ca*++Mg** contents of the 
water in the wet and dry seasons were greater. In the (Ca**+Mg**)/HCO, diagram, most of 
surface water and groundwater samples fell to the upper left of the 1:1 line, which indicated that 
the dissolution of carbonate played a dominant role in the formation of hydrochemistry. 
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Fig. 4 Ion ratio of different types of water body. (a), (Ca?*+Mg?*)/(HCO;+S0; ); (b), (Ca2*+Mg?")/SO; ; (c), 
(Ca?*+Me**V/HCO;. 


3.4 Distribution of heavy metals 


Cu, Mn, Al, Ni, Pb, Cu, and Ge were detected in groundwater samples retrieved from the HRB. 
On average, Al was the most concentrated heavy metal in groundwater samples, followed by Zn, 
Cu, Mn, Ni, Ge, and Pb. Analysis revealed 2.50, 1.73, and 2.45 times higher concentrations of Cu, 
Mn, and Al, respectively, in the mining areas and downstream of the HRB, relative to their 
concentrations in the unmined area. We produced a composition diagram for heavy metals 
detected in groundwater based on the Class III water standard stipulated in Groundwater Quality 
Standards (GB/T14848—2017, 2017) (Fig. 5). All groundwater samples collected from the study 
area were clearly within the upper limits defined in the Class III water standard, which suggested 
that groundwater in the study area posed a low risk to human health, and that domestic and 
industrial (e.g., mining) activity had little impact on the heavy metal contents of groundwater. 

3.5 Isotopic composition characteristics 

3.5.1 Precipitation 

5D and 8!8O values of water were relatively stable during the water cycle, with precipitation 
being a crucial step, and their variation in precipitation informed the evolution of the hydrological 
cycle. By calculating values of ôD and 8!°O in meteoric waters across the globe, Craig produced a 
global meteoric water line (GMWL) equation: 5D=86!8O+10. 5'8O and ôD analyses of the 


meteoric water in the HRB (Fig. 6) yielded the following local meteoric water line (LMWL) 
equation: 5D=7.605'80+4.41. As shown in Table 2, the values of 5D and 5'8O ranged from 
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—13.6%0 to —-6.67%o and from —101.44%p0 to —-45.52%o, respectively. Both the slope and intercept of 
the LMWL derived from the measurements were smaller than those of the GMWL, which was a 
result of the evaporation of migrating water vapour that encountered the dry ambient environment 
in the northwestern interior of China, where the HRB is located. Values of 5D and 6!8O in 
precipitation varied with time to some extent. Specifically, these two isotopes were enriched 
during the wet season and relatively depleted during the normal water and dry seasons. 
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Fig. 5 Composition of heavy metals in the groundwater. (a), Al; (b), Mn; (c), Ni; (d), Cu; (e), Zn; (£), Ge; (g), Pb. 
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Fig. 6 Relationship between ôD and 5!8O in precipitation. LWML, local meteoric water line; GWML, global 
meteoric water line. 


Table 2 5D and 8!°O data of precipitation 


5'0 (%o) 5D(%o) 
Season 
Maximum Minimum Average Maximum Minimum Average 
Normal water —8.46 —13.42 —10.67 48.11 -75.89 -59.85 
Wet —6.67 -10.94 -8.71 —62.42 -97.07 -78.05 


Dry -9.77 —13.60 -11.93 -73.60 —101.44 —86.70 
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3.5.2 Surface water, groundwater, and mine water 


Considerable variation was observed in 6'8O values of surface water and groundwater in the study 
area during the normal water, wet, and dry seasons, with larger variations found in surface water 
than in groundwater. 5!8O values of surface water ranged from —6.21%o0 to —8.5%o (with an 
average of —7.45%bo), from —4.33%o0 to —9.23%bo (with an average of —7.71%o0), and from —6.81%0 to 
—9.01%0 (with an average of —7.77%o) during the normal water, wet, and dry seasons, respectively. 
5!8O values of groundwater varied from —11.02%o to —7.57%o (with an average of —8.91%o), from 
—10.66%o0 to —7.23%o (with an average of —8.72%o), and from —10.35% to —7.5%o (with an average 
of —8.86%o) during the normal water, wet, and dry seasons, respectively. 5'8O value of both 
surface water and groundwater peaked during the wet season. 5'8O values of surface water and 
groundwater at most sites fell within the range of 5!°O in precipitation. This finding showed that 
atmospheric precipitation was the main source of surface water and groundwater. 

An analysis of the relationship between 5D and 5!°O in the unmined area (Fig. 7a) revealed a 
weak hydraulic connection between its surface water and groundwater. At most sites, 5D and &!80 
values were higher in surface water than in groundwater, which can be ascribed to pronounced 
isotopic enrichment of surface water due to heavy evaporation. The similar 6'8O and ôD values 
found in surface water and groundwater during the wet season suggested a relatively strong 
hydraulic connection between these values during this season. Figure 7b showed the relationship 
between 6'8O and 8D in different types of water body in the mining areas and downstream of the 
HRB. Compared with the unmined area, surface water and groundwater in the mining areas and 
downstream of the TRB were similar in isotopic composition, which suggested a stronger 
hydraulic connection, which may be attributed to an accelerated exchange between these types of 
water caused by human activity and coal mining. Even greater similarity in isotopic composition 
can be observed in surface water and groundwater during the wet season, reflected that frequent 
exchange between these types of water. Groundwater data for the dry season were closer to those 
from the LMWL, which was a consequence of intensified evaporation of shallow groundwater 
due to decreased precipitation. 5!8O and ôD values of the mine water did not significantly varied 
with sampling time. 5!8O and ôD values of the mine water were more consistent with those in the 
groundwater than with those of surface water. This finding suggested that the aquifer in the study 
area was deep, and had a weak hydraulic exchange capacity, and groundwater was a more capable 
recharge source for the mine water. 
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Fig. 7 Relationship between 5D and 8!8O. (al and a2), normal water season; (b1 and b2), wet season; (cl and 
c2), dry season. 
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Table 3 8D and 8!8O of different types of water body in the natural watershed 


Surface water 


Season 3580 (%0) 3D (%o) Fitting equation 
Maximum Minimum Average Maximum Minimum Average 

Normal water -6.21 -7.36 —6.97 -50.98 —61.79 -56.51 8D=7.875!8O-1.67 
Wet -6.61 —8.16 -7.39 -56.71 —64.38 —60.24 ôD=4.718'80-25.45 
Dry —6.81 -7.86 -7.23 -57.75 -59.74 -58.54 8D=1.135'8O-50.38 

Groundwater 

Normal water -7.57 —11.02 -9.21 -57.55 -86.19 —69.44 5D=8.285'80+6.85 
Wet -7.99 —10.66 -9.06 -56.71 —86.08 —69.42 öD=10.365!'80+24.37 
Dry -7.90 -9.98 -7.23 —62.31 -76.60 —67.12 öD=6.418'80-10.29 


Table 4 8D and 5!8O of different types of water body in the mining areas and downstream of the Hailiutu River 
Basin 


Surface water 


Season 35180 (%0) 3D (%o) Fitting equation 
Maximum Minimum Average Maximum Minimum Average 
Normal water -7.47 —8.50 -7.98 —58.40 —68.90 64.47 38D=4.315'8O0-30.02 
Wet 4.31 -9.23 -8.01 -53.51 -70.01 —62.57 ôD=2.538!'80—42.29 
Dry -7.28 -9.01 —8.30 -56.89 -73.30 —64.34 8D=6.255'8O-12.48 
Groundwater 
Normal water -7.61 —9.80 —8.60 -58.53 -77.46 —65.01 8D=7.385'8O-1.60 
Wet -7.23 -10.25 —8.37 -56.52 -80.95 —63.76 &D=6.805!8O-7.35 
Dry —8.05 —10.35 —8.85 -59.75 -78.58 —65.31 öD=7.218'80-1.51 
Mine water 
Normal water —8.85 -9.96 -9.30 —69.14 -81.86 -74.91 3D=8.605'80+4.80 
Wet —8.24 —10.20 -9.09 —63.28 -73.30 —67.37 5D=4.708'8O-24.70 
Dry -9.19 —10.84 -9.90 -77.11 -85.03 -80.91 5D=4.805'8O-32.90 


3.6 Transformation among precipitation, surface water, groundwater, and mine water 
3.6.1 Transformation of different types of water body in the natural watershed 


We analyzed transformation between surface water and groundwater based on the changes in TDS 
and 6!8O along the course of the HRB. Similar trends were observed in TDS and 8!8O values in 
the surface water and groundwater during the normal water, wet, and dry seasons (Fig. 8a). 
Surface water in the HR was enriched in !8O, while surface water in the YR and groundwater 
beneath it showed similar !8O values. The enrichment of '8O can be attributed to the intense 
evaporation of surface water in the unmined area, which was located in the Maowusu Desert. 
5!8O values of the groundwater beneath the HR first decreased and then increased along its course, 
whereas TDS did not significantly vary, indicating that the groundwater received recharge from 
surface water enriched in !8O. Overall, 5'!8O and TDS variations in surface water in the YR and 
groundwater beneath it along its course were not significant. 5'8O values of surface water in the 
YR gradually decreased along its course, suggesting that the surface water received recharge from 
groundwater depleted in '8O. However, the gradual increase in 5!8O in groundwater beneath the 
YR along its course indicated that the groundwater similarly received recharge from !8O-enriched 
surface water. Evidently, surface water in the YR and the groundwater beneath it intensively 
recharge one another. Multivariate linear mixed model was used to calculate the ratios of 
transformation among precipitation, surface water, and groundwater during the normal water, wet, 
and dry seasons (Fig. 8c). On average, precipitation contributed 15.6%, 18.7%, and 13.2% of the 
recharge received by surface water and 7.3%, 9.8%, and 5.2% of the recharge received by 
groundwater during the normal water, wet, and dry seasons, respectively. In addition, surface 
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water was responsible on averages for 44.9%, 51.0%, and 24.6% of the recharge received by 
groundwater and received 28.3%, 26.6%, and 33.5% of its recharge from groundwater during the 
normal water, wet, and dry seasons, respectively. These data showed that the predominant mode 
of water transfer between the surface water and aquifers in the unmined area involved the 
recharge of groundwater by infiltrated surface water during the normal water and wet seasons, 
and the lateral recharge of surface water by groundwater during the dry season. 
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Fig. 8 (al and a2), 5'8O and TDS variations in the surface water and groundwater in the unmined area along the 
course of the Hailiutu River Basin; (b), schematic diagram of water transfer between the surface and aquifers in 
the unmined area; (c), proportion of transfer between surface water and groundwater in the unmined area. 


3.6.2 Transformation of different types of water body in the mining areas and downstream of 
the HRB 


Figure 9a showed the variations in 8'8O and TDS values of surface water and groundwater in the 
mining areas and downstream of the TRB. The variations in 6'8O values in the surface water and 
groundwater in and around the mining areas (K1—H13) were complex and displayed no distinct 
patterns. TDS content was high in surface water and varied over a broad range, and TDS content 
of surface water was higher than that of groundwater. The overall content of TDS in the 
groundwater was low but increased, which indicated that groundwater was recharged by water 
sources with a high TDS content. Therefore, we speculated that groundwater was recharged by 
surface water with a high TDS content. In the lower reaches of the HR (H14—-h11), surface water 
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Fig.9 (al and a2), 8!8O and TDS variation in the surface water, groundwater, and mine water in the mining area 
and downstream of the Hailiutu River Basin; (b), schematic diagram of water transfer between the surface and 
aquifers in the mining areas and downstream of the Hailiutu River Basin; (c), proportion of transfer between 
surface water and groundwater in the mining area and downstream of the Hailiutu River Basin. 


was high in TDS and showed higher 5'8O values than those of groundwater, in which 6'8O and 
TDS values gradually increased along the course of the river, suggesting that the mode of transfer 
in this area involved the recharge of groundwater by surface water. The mine water was low in 
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5!8O and high in TDS at all the sites, which can be attributed to its great depth and the low flow 
rate of the aquifer. As seen in Figure 8a, surface water in the Yingpanhao mining area (K1, H5, h2, 
and h3), the Dahaize mining area (K2, H7, and h5), and the Balasu coal mine area (K3, K4, h6, h7, 
H10, H12, and H13) were generally high in TDS, indicating that dewatering and drainage 
operations significantly impacted surface water. 

Multivariate linear mixed model was employed to calculate the ratios at which precipitation, 
surface water, and groundwater recharge each other. As shown in Figure 8c, precipitation was 
responsible for the averages of 20.7%, 26.6%, and 15.8% of the recharge received by surface 
water, and 12.4%, 19.7%, and 7.9% of the recharge received by groundwater during the normal 
water, wet, and dry seasons, respectively. Compared with the unmined area, in the mining areas 
and downstream of the HRB, precipitation greatly contributed to the recharge received by both 
surface water and groundwater, particularly during the wet season. Surface water contributed 
47.8%, 61.1%, and 39.9% of the recharge received by groundwater, and received 36.3%, 31.7%, 
and 38.7% of its recharge from surface water during the normal water, wet, and dry seasons, 
respectively. Compared with those in the unmined area, the corresponding recharge proportions 
for the mining areas and downstream of the HRB were also high. Specifically, precipitation 
contributed 2.6%-—7.9% and 2.7%—9.9% more to the recharge received by surface water and 
groundwater, respectively; surface water contributed 2.9%—15.3% more to the recharge received 
by the groundwater; and groundwater contributed 5.1%—8.04% more to the recharge received by 
surface water. These results showed that factors such as coal mining and other types of human 
activity accelerated the transformation among precipitation, surface water, and groundwater. Mine 
water was composed primarily of the discharge resulting from drainage operations that occur 
during coal seam mining. Calculations revealed that precipitation, surface water, and groundwater 
contributed 4.1%, 15.3%, and 69.4% of the recharge received by mine water during the wet 
season, respectively. The proportions were 5.2%, 17.9%, and 72.4% for the normal water season, 
and 2.9%, 10.2%, and 61.3% for the dry season, respectively. Mine water was transferred into 
surface water and groundwater at proportions of 10.9%—40.1% and 8.3%-—30.21%, respectively. 
Compared with the unmined area, in the mining areas and downstream of the TRB, precipitation, 
surface water, and groundwater recharge each other at a higher proportion. In the mining areas 
and downstream of the TRB, mine water primarily received recharge from groundwater, while 
dewatering and drainage operations implemented during coal mining mainly had a significant 
impact on surface water. These findings were consistent with those derived from earlier 
hydrochemical analysis. Surface water and precipitation have minor impacts on mine water. The 
predominant mode of transfer in the mining area and downstream of the TRB involved the 
recharge of groundwater by surface water, which can be attributed to dewatering and drainage 
operations that occur during coal mining period. These actions lowered groundwater level in the 
mining areas, and accelerated the convergence of groundwater from the aquifer overlying the 
coal-bearing formation in the goafs, resulting in an increase in the infiltration of surface water 
into the ground. 


4 Discussion 


Generally, due to the continuous influence of river water evaporation, TDS flowing to surface 
water along the river will increase (Li et al., 2021), while the larger TDS values are concentrated 
in the middle and upper reaches of the HRB, and the changes in TDS and the main ions observed 
in the lower reaches were smaller, which may be due to the impact of coal mining in the middle 
and upper reaches of the HRB. Coal washing and drainage increase the contents of Ca*+, Mg”, 
and SO; in surface water. Except for the mining areas, TDS of the surface water increases along 
the river, but the overall variation is small, which shows that the impact of coal mining on the 
basin is mainly concentrated in the mining areas and has little impact on the non-mining areas 
(Huang et al., 2018). TDS of groundwater in the HRB changes only slightly from upstream to 
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downstream, indicating that coal mining and human activities have little impact on groundwater. 

Precipitation in the HRB had 8'8O values ranging from —13.6%o to —6.67%o, and 5D values 
ranging from —101.44%0 to —45.52%o. The slope of the LMWL equation derived from the 
measurements was smaller than that of the GMWL, because the study area is located in the inland 
of Northwest China, and local precipitation is strongly affected by evaporation. The water vapour 
in the atmosphere is affected by evaporation during movement, resulting in enhanced 
fractionation and dilution of hydrogen and oxygen isotopes in the remaining water vapour (Zeng 
et al., 2021). The hydraulic connection between surface water and shallow groundwater varied 
between seasons. A relatively strong hydraulic connection was apparent between surface water 
and groundwater during the wet season because the increase in precipitation in the wet season 
enhanced the mutual recharge between surface water and groundwater (Sun et al., 2021). Surface 
water and groundwater were more similar in isotopic composition in the mining areas and 
downstream of the HRB than in the unmined areas, suggesting that mining strengthens the 
hydraulic connection between surface water and groundwater. Due to the goafs caused by coal 
mining, the speed of seepage of surface water is accelerated (Jhariya et al., 2016). The isotopic 
composition of mine water in the study area was more similar to that of groundwater, indicating 
that mining activity disturbs the aquifers overlying the coal-bearing formation, and accelerates the 
downward recharge of groundwater (Miao et al., 2021). Based on multivariate linear mixed model, 
the results from mine water show that water discharged as a result of dewatering and drainage 
operations during coal mining is transferred to the surface at a higher ratio; however, there is a 
closer hydraulic linkage between mine water and groundwater. This result is possibly explained 
by the water drained from coal mine firstly discharging into the river, mixing with river water, and 
then seeping into the ground. In previous studies, the factors affecting the transformation among 
precipitation, surface water, and groundwater are complex, which mainly include climate change, 
topography, and human activities (Yang et al., 2017; Bouzekri et al., 2020). In the HRB, the 
impact of coal mining on the transformation among precipitation, surface water, and groundwater 
is the most pronounced in the mining areas, while water transfers between these compartments in 
the unmined areas were affected primarily by changes in seasons and topographic factors. 


5 Conclusions 


In the HRB, rock weathering governs the hydrochemical types of groundwater and mine water, 
while that of surface water is controlled by evaporation. Human activity (e.g., coal mining and 
agricultural production) causes considerable changes in the hydrochemical types of surface water 
in and around the mining areas, and leads to the most significant increases in concentrations of 
Nat and SOF. 

Precipitation is the main supply source of water resources in the HRB. In the unmined areas of 
the HRB, the primary mode of water transfer involves the recharge of groundwater by surface 
water during the normal water and wet seasons, and the recharge of surface water by groundwater 
during the dry season. In the mining areas and downstream of the HRB, year-round recharge of 
groundwater by surface water is the primary mode of water transfer. In the mining areas, 
precipitation contributed to the recharge received by surface water and groundwater by an average 
increases of 2.6%-—7.9% and 2.7%-—9.9%, respectively. Groundwater in the Salawusu Formation 
constitutes up to 61.3%—72.4% of mine water. 

This study revealed the transformation among different types of water body in the basin under 
the influences of natural watershed and coal mining qualitatively and quantitatively. However, the 
impact of coal mining on the transformation among different types of water body in the study area 
need to be analyzed with the long-term data. 
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